DFT study of graphene antidot lattices: The roles of geometry relaxation and spin 



o 
o 

(N 

< 

00 



Joachim A. Furst,^^EI Thomas G. Pedersen,^ Mads Brandbyge,^ and Antti-Pekka Jauho^'^ 

'DTU Nanotech - Department of Micro- and Nanotechnology, 
Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark 
^Department of Physics and Nanotechnology, Aalborg University, DK-9220 Aalborg East, Denmark 
"^Department of Applied Physics, Helsinki University of Technology, P. O. Box 1100, FI-02015 HUT, Finland 

(Dated: April 8, 2009) 

Graphene sheets with regular perforations, dubbed as antidot lattices, have theoretically been pre- 
dicted to have a number of interesting properties. Their recent experimental realization with lattice 
constants below 100 nanometers stresses the urgency of a thorough understanding of their elec- 
tronic properties. In this work we perform calculations of the band structure for various hydrogen- 
passivated hole geometries using both spin-polarized density functional theory (DFT) and DFT 
based tight-binding (DFTB) and address the importance of relaxation of the structures using either 
method or a combination thereof. We find from DFT that all structures investigated have band gaps 
ranging from 0.2 eV to 1.5 eV. Band gap sizes and general trends are well captured by DFTB with 
band gaps agreeing within about 0.2 eV even for very small structures. A combination of the two 
methods is found to offer a good trade-off between computational cost and accuracy. Both methods 
predict non-degenerate midgap states for certain antidot hole symmetries. The inclusion of spin 
results in a spin-splitting of these states as well as magnetic moments obeying the Lieb theorem. 
The local spin texture of both magnetic and non-magnetic symmetries is addressed. 

PACS numbers: 73.63.Fg, 73.63.-b 
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I. INTRODUCTION 

Graphene, the single-atom thick two-dimensional 
sheet of carbon atoms, has stimulated considerable 
experimental^ and theoretical research^ as well as pro- 
posals for future nanodevices^. Various graphene-based 
applications have been realized in recent years^ and the 
relevance for application in devices is heavily increased by 
the rapidly improving ability to pattern monolayer films 
with e-beam lithographj*i where features on the ten-nm 
scale have been obtained^J^,. Moreover, recent advances 
in chemical vapor deposition of graphene (see e.g. Ref. 
0) are promising for fabrication of large area, high qual- 
ity devices. 

Yet another way of nanoengineering graphene consists 
of defining an antidot lattice on graphene by means of 
a regular array of nanoscale perforations. This theo- 
retical idea was introduced by Pedersen et aZ.^'^^, who 
showed using tight-binding calculations that antidot lat- 
tices change the electronic properties from semimetal- 
lic to semiconducting with a significant and controllable 
band gap. Such structures have recently been realized 
experimentally by Shen et alM and Eroms et alJ^ with 
lattice spacings down to 80 nm. Quantum dots and 
graphene ribbons have been demonstrated with dimen- 
sions of only a few nmH. Very recently, Girit et alM have 
studied the dynamics at the edges of a growing hole in 
real time using a transmission electron microscope. Both 
in the experiment and in Monte Carlo simulations they 
find the zig-zag edge formation to be the most stable 
structure. This is in agreement with the findings of Jia 
et al}^ who demonstrate a method to produce graphitic 
nanoribbon edges in a controlled manner via Joule heat- 
ing. This opens the possibility of making antidot lattices 



with a desired hole geometry. 



FIG. 1: (Color online) The unit cell of the {4, 2}o(left) and 
{6, 5}a (right) system. The hexagonally shaped unit cells are 
repeated in plane to form a honeycomb lattice of antidots. 
The carbon atoms (green) are hydrogen terminated (white) 
along the hole edges. 

While carbon nanotubes, graphene and recently 
graphene ribbons have been studied extensively using 
first principles methods, antidot lattices in graphene have 
mainly been treated with simpler models^iii^. The very 
recent work by Vanevic et ali^ uses a 7r-orbital tight- 
binding model to study antidot lattices with rather large 
lattice constants (these systems are more easily accessed 
experimentally but cannot be analyzed in terms of ab ini- 
tio methods). Their focus is on the possible occurrence 
of midgap states without introducing defects in the an- 
tidot lattice, as was the case in the original proposal by 
Pedersen et al. 

Many studies on magnetization have been reported for 
various graphene structure o^^dSii9i20,2i,22 _ r^^iQ origin of 

the magnetism can be understood based on the theory 
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by Liebii, and the subsequent related work by Inui et 
alM- on the properties of the bipartite lattice. Single 
vacancies and their spin properties have been studied 
by, e.g. Lehtinen et al}^ and Palacios et al}^; the lat- 
ter paper also investigates voids in both graphene and 
graphene ribbons in detail using a mean-field Hubbard- 
model. Magnetization has also been studied in carved 
slits2£, finite ribbons^^ and flake o^^'^^ as well as rings^i 
and notchea^^. Recently, DFT treatments of magnetic 
properties of nano-holes in graphene^^ and graphene 
films^^ have been published. 

The realized antidot lattices with hole sizes of several 
tens of nanometers and even larger lattice spacings in- 
volve several thousands of atoms in a unit cell and are 
computationally too costly to be treated with DFT in 
a systematic manner. The DFT based tight-binding 
method, DFTB^*, however, allows one to address such 
large systems. The difference in computational cost be- 
tween DFT and DFTB is for the present study found to 
be at least a factor of thirty. We thus investigate the 
accuracy of DFTB compared to DFT on much smaller 
antidot lattices in terms of the band structures'^^. Since 
geometry relaxation is the most costly task in DFT we 
also investigate the cost benefits of combining the two 
methods. By using DFT and elaborating on the role of 
spin, we also wish to address some of the main features 
found specifically for antidot lattices on a tight-binding 
level of the theory. 

This paper is organized as follows. In Sec. [IT] we intro- 
duce the antidot lattice systems and the methods used. 
The equilibrium geometries and band structures obtained 
using both DFT and DFTB and a combination thereof is 
given in Sec. IIIII together with a detailed investigation of 
the spin properties. We conclude with a short summary. 



II. SYSTEMS AND METHOD 

The specific realization of the antidot lattice we con- 
sider in this paper is a hexagonal (triangular) array of 
holes in a graphene sheet as proposed by Pedersen et al.^ . 
Within the hexagonal unit cells there can be different hole 
geometries, and two examples of high symmetry holes 
are shown in Fig. [TJ Below, geometries are fully relaxed 
but as a starting point ideal geometries using fixed bond 
lengths and angles of 120° are constructed. These ge- 
ometries furthermore provide a straightforward notation 
for the structures. Thus, we designate antidot lattices 
with circular holes according to the notation {L,R}q, 
where L is the side length of the unit cell and R the 
hole radius, both measured in units of the graphene lat- 
tice constant a=2.46 A giving a C-C bond length of 1.42 
A. Similarly, for triangular holes, we apply the notation 
{L, D}/\, where D denotes the side length of the holei^. 
The holes are passivated with H using a C-H bond length 
of l.I A and consist almost entirely of zig-zag edges. 
These structures are idealized but may well be within 
experimental reach given the recent advancements^i^. 
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TABLE I: The band gaps for various systems calculated with 
either DFT or DFTB using geometries obtained with different 
methods for relaxation. All values are in eV. 



For the first principles calculations we have used the 
ab initio pseudopotcntial DFT as implemented in the 
Siesta code^^ to obtain the electronic structure and 
relaxed atomic positions from spin-polarized DFT—. 
We employ the GGA PBE functional for exchange- 
correlation'^". 

For the DFTB results^, we use the original (C,H) 
parametrization of Porezag et al^ which does not in- 
clude spin. 



III. RESULTS 

For all considered structures a structural relaxation 
with DFT leads to a shrinking of the hole, of the order 
of 1 %, resulting in C-C bonds close to the edges stretch- 
ing and contracting in the range of 1.39 - 1.45 (1.43 for 
{5, 3.5}q) a. a few bond lengths away from the hole 
edge the C-C bond length remains unaltered at 1.42 A. 

In the case of relaxation with DFTB the picture is quite 
similar. Edge-atom C-C bond lengths vary from 1.39- 
1.42 A for all systems but {6, 5}a, which has variations 
of 1.38 - 1.44 A. The shrinking of the hole size is smaller 
than 1%. 

The results for the band gaps for five different systems 
are summarized in Table IH The band structures are cal- 
culated using both DFT and DFTB on structures relaxed 
at different accuracy levels and thus at different compu- 
tational costs. The combinations are not exhaustive but 
represent a relevant set aimed at saving computational 
costs. The relaxation type is given in the second row of 
Table n DFTB is expected to match DFT better as L/R 
increases due to decreased importance of edge details. 
The systems chosen here are mostly edge-dominated (low 
L/R ratio) and thus represent the worst case scenario. 

Using DFT we find band gaps ranging from 0.2 to 
1.5 eV confirming that the antidot lattice turns the 
semimetallic graphene into a semiconductor-. However, 
only spin-polarized DFT predicts a band gap for the 
{6, 5}a structure which will be discussed in detail be- 
low. 

Pedersen et al^ demonstrated a scaling-law between 
the hole size and the band gap for large L/R ratios but no 
such simple picture for small L/R ratios emerged. This 
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FIG. 2: (Color online) The band gaps for DFTB, DFT on 
DFTB-relaxed geometry and DFT on unrelaxed geometry 
plotted versus pure DFT results. Points above (below) the 
dotted line are thus overestimated (underestimated) com- 
pared to pure DFT. Note, that DFTB calculates the elec- 
tronic structure without spin and fails to predict a band gap 
for the {6, 5}a structure. 



trend agrees well with our edge-dominated systems were 
no simple scaling between the hole size and the band gap 
could be identified. 

As illustrated on Fig. ^ DFTB in general gives a larger 
band gap than DFT (the circles lie consistently above the 
dashed line). This tendency is enhanced the more edge- 
like the structure becomes. On average, the DFTB gap 
is 20% larger than the DFT value for the four structures 
with circular perforations. The discrepancy increases for 
structures with holes occupying a large portion of the 
unit cell, such as {5,3.5}q. Moreover, for {6, 5}a the 
omission of spin effects leads incorrectly to a vanishing 
DFTB band gap in agreement with non-spin DFT. In the 
cases of the {4, 2}q and {6, 5}a systems the band struc- 
tures are shown in Fig. [3] left and right, respectively, cal- 
culated using DFT (DFTB) in the upper(lower) panel. 
We see that the shape of the bands corresponds qualita- 
tively for the two methods. 



A. Relaxation 

We next analyze the importance of relaxation. The 
clear trend is that relaxation increases the band gap. 
This is illustrated for DFT in Fig. [5] (crosses correspond- 
ing to the unrelaxed structure lie below the dotted line), 
as well as in Table |T] for DFTB. Comparing unrelaxed 
results with fully relaxed results we see from Table |T] a 
change in band gaps within 10 % and 15 % using DFTB 
and DFT, respectively. Only in the case of DFT does the 
effect of relaxation increase the more edge-dominated the 



system becomes. The change using DFT for the {4, 2}q 
system is 8 % compared to 16 % for {5, 3.5}q. 

It must be emphasized that larger differences between 
initial and relaxed geometries may well give rise to a 
larger discrepancy between their band gaps. However, 
even for the case of a single passivated edge-defect the 
difference in relaxed and unrelaxed DFT band gaps is 
less than 10 %. 



{4,2} 
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FIG. 3: (Color online) Band structures for the {4, 2}o (left 
column) and {6, 5}a (right column) systems using DFTB (up- 
per panels) and DFT (lower panels). The dotted blue curves 
in the lower panels are DFT results using structures relaxed 
with DFTB. The arrow pointing up (down) indicates that 
bands are filled with majority (minority) electrons only. 

By relaxing the geometry with DFTB the DFT re- 
sults are improved from the non-relaxed case as shown in 
Fig. 2 (the squares are closer to the dashed line than the 
crosses). Compared to pure DFT the largest difference 
in the band gap is again found for the edge-dominated 
{5, 3.5}q system: it is now 11 % compared to 16 % with- 
out relaxation. For the larger {6, 5}a structure we find 
the same values as for pure DFT. The DFT results are 
shown for both DFT-relaxed and DFTB-relaxed struc- 
tures on Fig. [3] indicated by thick and dotted line, re- 
spectively. The different geometries do not change the 
bands notably. 

Using DFT on DFTB-relaxed structures is thus an ap- 
proach with a good trade-off between accuracy and com- 
putational cost. This finding is of great practical use, 
since relaxation is very costly in DFT. 



B. Magnetic properties 

For the {6, 5}a system with both non-spin-DFT and 
DFTB there are three (one nearly doubly degenerate) 
bands with weak dispersion at zero energy. Introducing 



spin leads to a clear splitting of these bands, i.e., to the 
formation of a band gap. For a comparison of DFT and 
DFTB, see Fig. |3l right column. The three bands be- 
low (above) the Fermi level are half-filled by majority 
(minority) spin electrons and are thus completely spin- 
polarized. The size of the band gap is thus also an indi- 
cation of the robustness of the magnetic statc"^^ . 
The magnetic moments of the structures can be un- 




FIG. 4: (Color online) Amplitudes of certain important 
states. (a,b); degenerate midgap states, (c): non-degenerate 
midgap state of {6, 5} a , (d): highest filled bands of {6, 3.6}o • 
All states are calculated at the F-point. 

derstood as a consequence of graphene being a bipartite 
lattice in the nearest neighbor approximation as shown 
by Liebii. According to Lieb's theorem the total mag- 
netic moment can be written as M = Na — Nb where 
Na(b) is the number of atoms occupying the A{B) sites of 
the bipartite graphene lattice. Thus, if the angle between 
the zigzag edges is 0° or 60° the edge-atoms belong to the 
same sublattice, while they belong to different sublattices 
if the angle is 120° or 180°. Consequently, the hexago- 
nal hole is non-magnetic and the triangular is magnetic. 
This is consistent with a MuUiken analysis from the DFT 
calculations which shows a non-zero magnetic moment 
only for the {6, 5}a system of 3.00 per unit cell. By 
inspection of the geometry we indeed find Na — Nb—5- 
As a continuation of Lieb's work, Inui et ali^ showed 
that such sublattice imbalance results in there being 
Na — Nb midgap states with zero energy. We thus ex- 
pect a degeneracy of 3 of the low-dispersion bands in the 
{6, 5}a case. 

In Fig. [21 lower right panel, there is one largely dis- 
persionless band just below two almost completely de- 
generate bands which have some dispersion especially at 
the F-point. Such band curvature was also found for hy- 
drogenated graphene ribbons by Kusakabe et aL?-l The 




FIG. 5: (Color online) The diS'erence in majority spin and 
minority spin densities for the {6, 5}a structure. Blue (white) 
indicates surplus of majority(minority) spin. 



F-point states for each band are shown in Fig. 31 where 
the strong localization is seen for the lowest band. Fig. 
mjc), whereas the bands with curvature, Fig. |llja,b), yield 
less localized states. Note also the alternation in the am- 
plitudes of the states between sublattices as proposed 
by Inui et al^ The state of the highest occupied (spin- 
degenerate) band for the non-magnetic {6, 3.6}q struc- 
ture is shown for comparison in Fig. [DJd). The splitting 
as well as the curvature is less pronounced for the unfilled 
states above the Fermi level showing particle-hole asym- 
metry. This asymmetry is expected due to the breaking 
of the symmetry of the bipartite lattice partly due to 
the DFT treatment beyond nearest neighbor as well as 
the passivation of the edges which changes the on-site 
potential at edge-sites. This is inherently also the case 
for DFTB. By inspection of the SIESTA HamiltonianSa 
we find an increase in on-site energy for passivated edge 
atoms as compared to atoms far from the edge. Vanevic 
et ai— find that a potential shift on the edge-atoms 
mainly causes a lifting of the degeneracy of the fiat bands, 
consistent with our observations. 

As mentioned above, the global spin is given by the 
sublattice imbalance. This does not, however, determine 
the local spin. For the hexagonally shaped hole struc- 
tures we find not only a zero global spin, but also a 
zero spin on all atoms. This explains the identical band 
gaps found using DFT with and without spin. Such non- 
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magnetic solutions have been found for finite graphene 
ribbons or graphene flakes by Jiang et al^ using DFT. 
They find a sudden transition from non-magnetic to mag- 
netic solutions going from a system size of [3,3] and [4,3] 
with numbers indicating rings in the graphene lattice 
along zig-zag and armchair directions. Such transitions 
are also seen in slits cut in graphene in a study by Ku- 
mazaki et alP''^ Viewing each edge in our hexagonal struc- 
tures as ribbons, the largest ribbon corresponds to a [4,3] 
ribbon ({6, 3.6}q structure). We thus expect local mag- 
netization to appear for slightly larger systems. 

For {6, 5}a we have the strongest polarization at the 
middle of each edge with maximum magnetic moment per 
atom being 0.24 /is . We note that each corner atom has a 
magnetization -0.03 /is. The edge atom magnetization is 
below the maximum of 1 /3 for graphene ribbons when 
the width becomes too large for edge-edge interactions. 
Our systems thus have edge-edge interactions which is 
expected due to the ribbon width of 6 rings. A plot 
of the spin-polarized density^^ is shown in Fig. [5l The 
majority spins reside mostly on the edges of the dominant 
sublattice sites. Neighboring sites on the other sublattice 
have minority spin-polarization. Note also the non-zero 
spin of the atoms far from the edges indicating interaction 
between neighboring hole-edges. 



ation using either method is found to increase the band 
gap with maximally 15%. Combining the two methods 
by performing a DFT-calculation on a DFTB-relaxed 
structure is found to give a good trade-off between accu- 
racy and computational cost facilitating the treatment of 
larger systems. However, even for unrelaxed geometries 
we find qualitative agreement with the DFT-relaxed ge- 
ometries. Trends for ideal geometries as presented here 
can thus be investigated without any relaxation in the 
non-magnetic case. 

Certain geometries are shown with DFT to have a 
non-zero total magnetic moment which is understood via 
Lieb's theorem as a consequence of sublattice imbalance. 
For these structures a spin-polarized treatment is needed 
to achieve even qualitative results for the band gaps. Lo- 
cal spin, which can occur regardless of the total mag- 
netic moment, is not observed for very edge-dominated 
systems. Sublattice imbalance leads to the occurrence of 
low-dispersion midgap bands. We find a lifting of degen- 
eracy of these otherwise degenerate bands on a perfect 
bipartite lattice as well as a considerable spin-splitting. 
These completely spin-polarized states are primarily lo- 
cated at the hole-edges. 
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IV. CONCLUSION 

Using DFT and DFTB we have calculated band gaps in 
various antidot lattice geometries. The computed band 
gaps range from 0.2 to 1.5 eV. In general, DFTB gives 
larger band gaps than DFT with the largest difference 
for non-magnetic structures of 44 %. Geometry relax- 
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We note, that DFT is known to underestimate the band 
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The systems are relaxed using the conjugate gradient (CG) 
method with a force tolerance of 0.01 eV/A. All atoms 
are allowed to move during relaxation but the unit cell 
is kept fixed with a vacuum between graphene layers in 
neighboring cells of 10 A. The mesh cutoff value defining 
the real space grid used is 175 Ry and we employ a double- 
polarized (DZP) basis set. A Monkhorst-Pack grid of 
(2,2,1) was found sufEcient for all structures. 
We relax the two outermost rows of atoms along the 
perimeter of the hole until the total energy per atom was 
converged to better than 10^* cV. 

Due to the larger cell used in this calculation we use less ac- 
curate settings (single-^" basis, 100 Ry mesh, one fc-point). 
This changes the magnitude of the spin somewhat as com- 
pared to the smaller cell calculations. The maximum mag- 
netic moment per atom is 0.27 /Ub for the large cell versus 
0.24 /is for the smaller. The qualitative picture of the spin 
distribution is, however, the same for both calculations. 



